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ABSTRACT 

We report the detection of variable linear polarization from Sgr A* at a wavelength of 3.5 mm, the 
longest wavelength yet at which a detection has been made. The mean polarization is 2.1 ± 0.1% at 
a position angle of 16 ± 2° with rms scatters of 0.4% and 9° over the five epochs. We also detect 
polarization variability on a timescale of days. Combined with previous detections over the range 150- 
400 GHz (750-2000 ^m), the average polarization position angles are all found to be consistent with a 
rotation measure of —4.4 ± 0.3 x 10^ radm~^. This implies that the Faraday rotation occurs external 
to the polarized source at all wavelengths. This implies an accretion rate ~ 0.2 — 4 x 10^* Mq yr^^ for 
the accretion density profiles expected of ADAF, jet and CDAF models and assuming that the region 
at which electrons in the accretion flow become relativistic is within 10i?s- The inferred accretion 
rate is inconsistent with ADAF/Bondi accretion. The stability of the mean polarization position 
angle between disparate polarization observations over the frequency range limits fluctuations in the 
accretion rate to less than 5%. The flat frequency dependence of the inter-day polarization position 
angle variations also makes them difficult to attribute to rotation measure fluctuations, and suggests 
that both the magnitude and position angle variations are intrinsic to the emission. 
Subject headings: galaxies: active — Galaxy: center — polarization — radiation mechanisms: non- 
thermal 



1. INTRODUCTION 

Linear polarization can be an important diagnostic 
of relativistic jets and accretion flows associated with 
black hole systems. In the case of the massive black 
hole in the Galactic Center, Sgr A*, the properties of 
its mm-wavelength linear polarization probes the accre- 
tion environment on scales inaccessible with other tech- 
niques (Bower et al. 1999a, b; Aitken et al. 2001; Bower 
et al.2001; Bower et al.2003; Marrone et al. 2006). The 
apparent absence of linear polarization at wavelengths 
exceeding 2.7 mm and sharp rise in polarization fraction 
at shorter wavelengths sets an upper limit to the rota- 
tion measure (RM). This limits the mass accretion rate to 
- 1O~'^M0 yr"i at distances of 10 - 1000 Schwarzschild 
radii from the black hole, which eliminates certain classes 
of accretion flow (Quatacrt & Gruzinov 2000; Agol 2000), 
but is consistent with CDAF and jet models (e.g. Falcke, 
Mannheim and Biermann 1993). The RM measures the 
accretion rate by serving as a proxy for the electron col- 
umn density once coupled with assumptions about the 
magnetic field. Equipartition between kinetic, magnetic 
and gravitational energy is often assumed to relate the 
electron density and magnetic field (e.g. Bower et al. 
1999a; Meha & Falcke 2001; Marrone et al. 2006). 

The discovery of variations in both the polarization an- 
gle (Bower et al. 2005) and fraction (Marrone et al. 2006) 
suggests two potential sources of variability. The varia- 
tions may be intrinsic and, therefore, offer evidence on 
the nature and structure of the emission region on a scale 
of ^ 10i?s. Changes in the RM along the line of sight 
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may also induce external polarization variability. Struc- 
ture in the accretion region on all scales within the Bondi 
radius can contribute to RM variations. There is thus 
strong motivation to accurately characterize the RM, its 
fluctuations, and the intrinsic variability of the polarized 
source. 

In §2 we present the detection of linear polarization 
of Sgr A* at 3.5 mm, showing that it varies on short 
timescales as well as relative to historical non-detections. 
Combining this in §3 with other data, our detection 
yields the best constraints so far on the RM. We dis- 
cuss the nature of the accretion flow and on the intrinsic 
source properties in §4. 

2. OBSERVATIONS 

Sgr A* was observed with the BIMA array at 3.5 mm 
between 28 March and 01 April 2004 UT. All five obser- 
vations were obtained in identical LST time ranges that 
corresponded roughly to 11 to 16 UT. Observations were 
obtained in a polarimetric mode that produced Stokes /, 
Q, U and V images, with the two frequency bands cen- 
tered on 82.76 and 86.31 GHz, each with 800 MHz band- 
width. Data were calibrated assuming that the flux den- 
sity of J 1733- 130 was steady at its mean value of 2.7 Jy. 
J1733-130 was observed for 5 minutes every 30 minutes. 
Observations were obtained in the B configuration giving 
a resolution of 8.6 x 2.6" in position angle (p. a.) 7°. The 
data were filtered to remove spacings shorter than 20 kA, 
similar to the processing of previous Sgr A* polariza- 
tion experiments. This u-v distance cutoff gives a total 
flux estimate accurate to better than 10% from 8.4 to 86 
GHz (Bower et al.2001). A priori amplitude calibration 
using system temperature and default gain information 
was applied to J1733-130. Both sources were phase self- 
calibrated. Polarization leakage terms from observations 
of 3C279 on 11 November 2003 were applied. 

Numerous sources of error contribute to the accuracy 
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with which we can measure the polarization fraction and 
p. a.. In addition to statistical errors that include noise 
from the sky and the telescope, the leakage of polar- 
ization from one polarization handedness to the other 
is an important effect. This converts a fraction of the 
total intensity into a false polarization signal. Leak- 
age calibration uses observations of a bright calibrator 
source and simultaneously solves for the polarization of 
the source and the leakage terms appropriate for indi- 
vidual antennas. Imperfect correction leaves a residual 
false polarization signal and introduces errors that exist 
in the calibrator observation. BIMA leakage solutions at 
1.3 mm were found to be stable over periods of months 
to years; variability in the solutions were due to changes 
in receiver orientation (Bower, Wright & Forster, 2002). 
Typical variations in antenna leakage terms were ~ 1.5% 
over 2 years, indicating a ~ 0.5% variation in the aver- 
age leakage solution for the 9 BIMA antennas. Solutions 
at 3.5 mm are expected to have the same characteristics 
although this has not been studied in depth. In addi- 
tion to variations in the leakage error, statistical errors 
in the calibrator data introduce errors. A comparison of 
actual measured p.a.s of 3C 279 at 3.5 mm and 1.3 mm 
from BIMA with measured p.a.s at 1.3 cm and 0.7 cm 
from the VLA indicate excursions of 10°- 20° (Bower et 
al. 2002). These deviations exceed those expected from 
leakage corrections and statistical errors. They are pos- 
sibly real source effects but may also be uncorrected in- 
strumental errors. 

The average flux densities of Sgr A* and linear polar- 
ization E- vectors are listed in Tabled In the polarization 
fraction and p. a. errors we include the effect of a 0.5% 
polarization leakage error. This error propagation is ac- 
curate in the case of the high SNR detections in the LSB 
but presents a lower limit to the errors in the low SNR 
detections in the USB. Intcr-day flux density variations 
are detected at a significant level, with a reduced-x^ of 
4.9 for the hypothesis of constant flux density for Sgr 
A* for 4 degrees of freedom. The mean flux density is 
1.93±0.07 Jy and the rms variability amplitude is 0.18 Jy. 

Polarization was detected in all five epochs, making 
this the first detection of polarization at 3.5 mm, the 
longest wavelength at which polarization has been de- 
tected in Sgr A*. The mean polarization fraction is 
2.1 ± 0.1% and its mean p. a. is 16 ± 2°. The polarization 
fraction and p. a. vary from day to day. The detection of 
polarization is more robust in the LSB (82.8 GHz) than 
in the USB (86.3 GHz). The origin of this difference is 
unclear. It may be due to changes in the leakage correc- 
tions and/or poor phase stability at the higher frequency. 
We have previously seen data from BIMA in which polar- 
ization results in one sideband were more reliable than 
in the other (Bower et al. 2001). The bandwidth used 
in these observations is sufficiently small that the linear 
polarization is not depolarized by rotation of the polar- 
ization vector across the observing band, as discussed 
below. 

3. THE ROTATION MEASURE 

Table |21 lists the mean polarization p.a.s for all pre- 
vious detections of linear polarization in Sgr A*. Due 
to systematic errors from polarization leakage calibra- 
tion, the poor quality of upper sideband data at 3.5 mm, 
and the relatively small spacing between the two side- 
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The daily average flux densities of Sgr A*, /err is 

DETERMINED FROM THE FLUX DENSITY SCATTER DETERMINED ON 
SHORT TIMESCALES FOR EACH DAY. cr IS THE STATISTICAL ERROR 
FOR Q AND U ; ERRORS IN THE POLARIZATION FRACTION AND 
P. A. INCLUDE THE EFFECTS OF 0.5% LEAKAGE ERROR. 
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TABLE 2 

Mean polarization fractions and p.a.s from measurements 

OF "AlTKEN ET AL. 2001,''BOWER ET AL. 2003 & 2005, AND 

'^Marrone et AL. 2006. The quantities ap and cr^ denote 

THE RMS variation IN THE POLARIZATION FRACTIONS AND P.A.S 

RESPECTIVELY EXCEPT IN THE CASE OF THE SINGLE-EPOCH 
MEASUREMENTS, WHERE THEY DENOTE THE ESTIMATED ERROR OF 
THE MEASUREMENT. 



bands, we are unable to compute a meaningful RM es- 
timate from the 3.5 mm data alone. The combined av- 
erage polarization data are consistent with a single RM 
of —4.4 ± 0.3 X lO^radm"^ and an intrinsic polariza- 
tion p.a. xo = 168 ± 8°. Fig. 1 plots the best RM fit 
to the mean polarization p.a. data. The best fit, with a 
rcduccd-x^ of 29 computed using the errors quoted in Ta- 
blelU is obtained if the 85 GHz p.a. is rotated by —180°. 
A fit to the unwrapped 85 GHz p.a. yields a reduced-x^ 
of 54 and RM = -2.0 ±0.4 x lO^radm"^, while an extra 
-180° wrap yields RM = -6.9 ±0.4 x lO^radm'^ with 
an associated reduced-x^ of 64. 

The reduced-x^ for the best fit is unacceptably high be- 
cause of the unrealistically small errors associated with 
the Aitken et al. (2001) observations. These single-epoch 
low-resolution measruements predate the identification 
of polarization variability in Sgr A*. We therefore fit 
only to interferomctric observations spanning multiple 
epochs. The resulting fit properties depend on whether 
one regards measurements from the two sidebands at 
v > 216 GHz, where available, as being mutually inde- 
pendent. If so, the reduced-x^ for a fit including the two 
-180°-wrapped 3.5 mm points is 0.9, but is 3.6 for the 
unwrapped 3.5 mm polarization pas. If not, the two x^ 
values are 0.8 and 2.1 respectively. In both cases, the so- 
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lution involving the wrapped 3.5 mm points is preferred. 
We regard these cases as bounding the true significance 
of the RM. For the case of independent sideband mea- 
surements, the associated best fit is xo = 163 ± 2° and 
RM = -4.38 ± 0.06 x 10^ radm^^. 

We adopt the RM = -4.4 ±0.3 x lO^radm"^ derived 
from all of the multi-epoch data. Given the variability of 
the p. a. and the uncertain systematics of different mea- 
surements, we consider the use of all data to provide the 
most conservative estimate of the RM. 
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Fig. 1. — The mean polarization p. a. of Sgr A* as a function 
of wavelength. The 85 GHz points have been derotated by —180°. 
Diamonds denote measurements for which polarization variability 
is detected, and their error bars denote the standard error of the 
mean of the p. a. variations. Triangles denote single epoch mea- 
surements only. The dashed line shows the best fit to all data, the 
solid line to all points excluding the Aitken et al. (2001) data, and 
the dot-dashed line to the unwrapped 85 GHz points, denoted by 
stars. 
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Fig. 2. — The accretion rate implied by our measurement of the 
RM for various accretion models using A/(,ij = 2.6 X 10^ Mq (Ghez 
et al. 1998). Two choices for rout are shown: one in which the outer 
scale is large (effectively infinite) and another in which it is only 
three times larger than the inner cutoff radius. 



4. DISCUSSION 

The accretion rate implied by this RM depends on den- 
sity and magnetic profiles assumed. Following the pre- 
scription outlined by Marrone et al. (2006),'* we write 
n oc r"^, for rin < r < rout and assume equiparti- 
tion between magnetic, kinetic and gravitational energy. 
The radius rin designates the point at which the flow 
becomes relativistic, and so ceases to contribute to the 

Note the typographical error in eq. (9) of Marrone et al. (2006) 

— 7/4 

in which one should have RMoc rj^ 



RM because its contribution to Faraday rotation is de- 
creased by a factor 7^/ log 7. The outer scale is a proxy 
for the coherence length of the magnetic field fluctua- 
tions. The estimated accretion rate for various ADAF 
{13 = 3/2) and CDAF (/? = 1/2) models is shown in Fig- 
ure |21 as a function of Tin. The accretion rate is in the 
range 0.2 - 4 x IQ-^Moyr-i for 2 Rs < < 10i?s, 
rout > 3nn and 1/2 < ^ < 3/2. 

The fact that a single RM accounts for the frequency 
dependence of all mean polarization p.a.s implies that the 
Faraday rotation occurs external to the polarized source. 
Internal rotation would cause the RM to vary as a func- 
tion of frequency. The emission is optically thick at all 
frequencies at which polarization is detected (Falcke et 
al. 1998; Zhao et al.2003, Yuan et al.2003), so if any 
internal Faraday rotation did occur, the diminuation of 
opacity effects with frequency, which increases the depth 
down to which one observes emission, would cause a cor- 
responding increase in the Faraday rotation path length. 

An interpretation of the p. a. jitter observed at 85 GHz 
in terms of RM variability would imply rms deviations of 
1.2 X lO^'radm^^, far smaller than the ~ 2 x lO^radm"^ 
rms deviations implied by the 340 GHz Marrone et 
al. (2006) fluctuations interpreted similarly. However, 
the absence of a clear frequency dependence in the inter- 
day jitter makes it difhcult to ascribe to RM fluctuations, 
and hence variations in the accretion rate. Table 2 shows 
the rms p. a. deviations for the multi-epoch observations 
at 85, 216, 230 and 340 GHz. These are inconsistent with 
the ly^^ dependence expected of RM fluctuations from a 
magnetoionic medium external to the source, suggesting 
instead that the jitter reflects changes in the intrinsic 
source polarization p. a.. Nonetheless, the fact that the 
fluctuations at these frequencies were not observed si- 
multaneously, coupled with the short time span of the 
85 GHz observations compared to the > 2 month - albeit 
sporadic - sampling of the 230 and 340 GHz measure- 
ments, still admits the possibility that the p. a. disper- 
sion observed at 85 GHz is unrepresentative of its long- 
term average. This appears unlikely. If the ~ 20° rms 
p. a. variations observed at 216-230 GHz were associated 
with RM fluctuations we would expect w 150° fluctua- 
tions at 85 GHz and would not expect to find x ^ over 
a set of disparate measurements. Moreover, the presence 
of intrinsic p. a. changes is unsurprising given that the 
polarization fraction is also intrinsically variable, as dis- 
cussed below. 

The lack of p. a. jitter attributable to RM fluctuations 
can be interpreted in terms of an upper limit in accre- 
tion rate fluctuations. The absence of clearly identifiable 
inter-day RM fluctuations suggests it is uniform on inter- 
day timescales. The consistency of the observations over 
a large number of disjoint epochs and frequencies with 
a single RM further suggests that the underlying accre- 
tion rate is constant on the timescale over which the ob- 
served polarization pas were averaged. The uncertainty 
in our RM fit places an upper bound on the accretion 
rate variations using RM oc M'^/^, valid under the as- 
sumption of equipartition between magnetic, kinetic and 
gravitational energy (e.g. Marrone et al. 2006). The 7% 

uncertainty in the RM implies Al fluctuations less than 
5%. The limit on M variations is largely consistent with 
the limits imposed by source flux density variations. In 
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the jet model the flux density scales as M^^^^^ (Falcke 
et al. 1993). The standard deviation of the 3.5 mm fluxes 
is 10%, comparable to the errors on the individual mea- 
surements, which imposes an upper limit of 14% on M 
fluctuations. Marrone et al. (2006) detect 10% rms inten- 
sity variations at 340 GHz. The five intensity measure- 
ments from Bower et al. (2005) at 216-230 GHz exhibit 

29% modulations, implying M fluctuations of 40%. Vari- 
ability in polarization fraction is detected by all multi- 
epoch observations (Fig. 01, at 85, 230 and 340 GHz, and 
presumably occurs at 112 GHz, where it was not detected 
at the 1.8% (l-a) level by a previous search (Bower et 
al. 2001). A previous limit of 1% linear polarization 
at 86 GHz (Bower et al. 1999b) demonstrates that it also 
varies on long timescales. 

Since the Faraday screen is external to the source the 
polarization amplitude fluctuations must be intrinsic. In- 
strumental depolarization effects are too low to explain 
the variability: bandwidth depolarization is only impor- 
tant for RMs greater than 2.7 x 10''radm~2 at 85 GHz, 
while beam depolarization is similarly improbable given 
the sub-mas size of Sgr A* at mm wavelengths (Bower 
et al. 2004; Shen et al. 2005). Spatial variation in the 
RM across the transverse extent of the source could de- 
polarize the emission, but at even 85 GHz this requires 
RM fluctuations ^RM > 7 x lO^radm'^ (Quataert & 
Gruzinov 2000). 

Both the high variability of the mm and sub-mm emis- 
sion (Zhao et al. 2004; Wright & Backer 1993; Tsuboi et 
al. 1999) and the linearly polarized emission of Sgr A* 
are possibly associated with its excess sub-mm emission 
(Serabyn et al. 1997; Falcke et al. 1998; Meha et al. 2001). 
ADAF models that fit the cm to-sub-mm spectrum in- 
clude at least two distinct populations of radiating par- 
ticles (Yuan et al. 2003), with the second important at 
V ^ 100 GHz in order to explain the sub-mm bump. 

The model of Yuan et al. (2003), in which the sub-mm 
emission is dominated by thermal electrons, overpredicts 
the polarization fraction at 85 GHz (cf. Fig.|2Jl. In this 
model the degree of linear polarization ranges from 32% 
at 85 GHz to 70% at 400 GHz assuming a uniform mag- 
netic field and that Faraday depolarization intrinsic to 
the source is unimportant (see their Fig. 3b). In the con- 
text of this model the ratio of the predicted to observed 
polarization levels can only be attributed to magnetic 
field inhomogeneity intrinsic to the source. This ratio is 
a factor of 3 higher at 85 GHz relative to the fraction in 
the range 150 — 400 GHz, over which it is constant within 
the errors. It is hard to account for such an increase 



in magnetic field inhomogeneity at 85 GHz, particularly 
if the source size only scales oc z^^^, as expected if its 
emission is self-absorbed. On the other hand, including 
synchrotron self-absorption effects, Goldston, Quataert 
& Igumenshchev (2005) show that the polarization frac- 
tion is expected to increase by a factor of three over the 
range 85 — 200 GHz. The quasi-spherical accretion polar- 
ization model of Melia, Liu & Coker (2000) and the two- 
component model of Agol (2000) predict a 90° p. a. flip at 
^ 280 GHz which is at variance with the measured p.a.s 
at 150, 230 and 340 GHz. 

We have reported here the detection of linear polariza- 
tion in Sgr A* at 3.5 mm. This enables us to calculate the 
rotation measure and set a limit on the accretion rate. 
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Fig. 3. — The mean polarization fraction of Sgr A*. The error 
bars plotted are those of Table 1^ Previous 86 and 112 GHz non- 
detections are marked with arrows. The error bars in the Aitken et 
al. (2001) measurements, marked with triangles, reflect uncertainty 
in the contribution from dust emission rather than variability as- 
sociated with the source. 



The lack of frequency dependence for position angle fluc- 
tuations indicates that they are intrinsic to the source. 
Our result favors RIAF/CDAF accretion models, with 
a shallow density distribution, over ADAF and Bondi- 
Hoyle accretion flows, which have a steep profile and are 
more likely to produce rapid RM variations. Future wide 
bandwidth simultaneous observations with CARMA and 
the SMA will fully characterize intrinsic and extrinsic 
changes in the polarization properties of Sgr A* and al- 
low us to investigate the accretion environments of other 
nearby low luminosity AGN, such as M81* (Brunthaler, 
Bower & Falcke 2006). 
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